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ABSTRACT 



We address the problem of identifying remnants of satel- 
lite galaxies in the halo of our galaxy with Gaia data. We 
make use of N-body simulations of dwarf galaxies be- 
ing disrupted in the halo of our galaxy combined with a 
Monte Carlo model of the Milky Way galaxy. The models 
are converted to a simulated Gaia catalogue containing a 
realistic number (~ 10^-10^) of stars. The simulated cat- 
alogue can be used to study how to handle the large data 
set that Gaia will provide and to study issues such as how 
to best retrieve information on substructure in the Galac- 
tic halo. The techniques described are applicable to any 
set of N-body simulations of (parts of) the Galaxy. 
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experience with analysing and visualising the enormous 
volume of information that will be present in the Gaia 
data-base. 

To achieve this we built a Monte Carlo model of the 
smooth components of the Galaxy containing a realistic 
number of stars and performed tree-code A/^-body simu- 
lations of satellites that are disrupted while orbiting the 
Galaxy. The Galaxy and satellite models were subse- 
quently combined and Gaia observations were simulated. 
In the following we summarise how the simulated Gaia 
survey was generated, with emphasis on the proper com- 
bination of the satellite and Milky Way models, and we 
present some exa mples of the simulated Gaia data. Fo r 
details we refer to lBrown. Velazquez. & AguilaJ (l2004l) . 



2. MONTE CARLO MODEL OF OUR GALAXY 



1. INTRODUCTION 



Cosmological models predict that halos of luminous 
galaxies are assembled through merging of smaller struc- 
tures. Remnants of past mergers will survive for a long 
time in the halo as debris streams because of the very long 
dynamical time scale. The Gaia mission offers a unique 
opportunity to search for and study these remnants with 
full phase space information in our galaxy's halo. How- 
ever disentangling the possibly many remnants from each 
other and from the background of Galactic stars will be 
very challenging. The Gaia catalogue will contain about 
one billion objects of which the stars in the individual 
debris streams may form only a very small fraction. In 
addition the streams are spread out all over the sky. Iden- 
tifying them will require a combination of search meth- 
ods that make use both of conserved dynamical quanti- 
ties of the debris stream (such as energy and angular mo- 
mentum, see e.g.jHelmi etIaIl|1999HHelmi & de Zeeu^ 
[2000) and astrophysical properties of the constituent stars 
(from photometric data). 

The goals of this work are: 1. study the retrieval from 
Gaia data of remnants of satellite galaxies that have been 
disrupted in the potential of our Galaxy; 2. include a real- 
istic model of the Galactic background population against 
which the remnants have to be detected; 3. gain practical 



The Milky Way model consists of three spatial compo- 
nents: a bulge with a Plummer density law, a double ex- 
ponential disk, and a flattened halo (c/a = 0.8) for which 
the density drops as r~^-^. The kinematics are modelled 
in a very simple manner: each component rotates with 
constant velocity dispersion, for the bulge an isotropic 
dispersion is assumed, for the disk a different velocity el- 
lipsoid is used for each spectral type (OBAFGKM), while 
for the halo the dispersions are the same for each spectral 
type. The individual stars are assigned an absolute mag- 
nitude My and a colour {V — I) from a Hess-diagram 
which is considered fixed for all Galactic components and 
all positions throughout the Galaxy. The Hess-diagram is 
taken from Table 4-7 of Mihalas & Binney (1981) and 
provides the relative numbers of stars in bins of abso- 
lute magnitude (My) and spectral type. These numbers 
integrated over spectral type provide the luminosity func- 
tion. Although this is a highly simplified and certainly not 
self-consistent model of the Galaxy it is good enough for 
providing the 'background' distribution in phase space 
against which the debris streams have to be found. 

For a magnitude limited survey, as is the case for Gaia, 
a straightforward Monte Carlo realisation of the Galac- 
tic model is potentially a very wasteful procedure. Most 
simulated stars will be too faint to be included in the sur- 
vey. Hence we used a strategy that minimizes wasted ef- 
fort by generating stars only within the magnitude limited 
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Figure 1. Intrinsic (circles, left vertical scale) and 
weighted ( triangles, right vertical scale) luminosity func- 
tions for our Milky Way model plotted vs. My. The 
weighted luminosity function is plotted in terms of the 
number of stars expected in the simulated Gaia survey. 



sphere centred on the observer. This results in a Galaxy 
model with a luminosity function that is weighted by the 
space density integrated over the volume limited sphere 
for each spectral type. This luminosity function is shown 
in Fig. [I] Nevertheless the generation of the Monte Carlo 
model is still a large computational task and in order to 
save time we decided not to simulate the part of the sky 
within Galactic coordinates: -90° < £ < +90° and 
—5° < b < +5°. Tracing the debris streams in this part 
of the sky will be difficult in practice due to the large ex- 
tinction in those directions. In our Galactic model ~ 80% 
of the stars lie in this region of the sky as seen from the 
Sun for a survey limited at F = 20. The resulting simu- 
lated survey of the Galaxy is fully sampled and contains 
3.2 X 10^ stars. For each of these we generated the 6 
phase space coordinates (positions and velocities), an ab- 
solute magnitude, a population type (bulge, disk or halo) 
and a spectral type. 



3. N-BODY MODELS OF DWARF GALAXIES 



The models of the debris streams were generated by 
simulating the disruption of dwarf galaxies orbiting our 
galaxy. The Milky Way is represented by a rigid poten- 
tial which is derived from a mass-model consisting of: a 
double exponential disc, a spherical bulge with Hernquist 
profile, and a logarithmic Halo potential with flattening 
c/a = 0.8. The dwa rf galaxies are represented by King 
models (j^ESIUll) with a mass of 5.6 x 10^ M© or 
2.8 X 10^ M0 and 10^ particles. The tidal radius and 
concentration parameter c = log (rt/rc) are 3150 pc and 
0.9, respectively. The satellites are placed on five differ- 
ent orbits which vary in apocentre, pericentre and the ini- 
tial inclination with respect to the Milky Way's disc. The 
simulations were evolved with a tree-code for ~ 10 Gyr. 
Multiple debris streams can be simulated by combining 
different A/'-body snapshots (different orbits and/or ages), 
and each snapshot can be rotated around the Z-axis or 
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Figure 2. Satellite orbits for models no. 1 and 4 pro- 
jected onto the principal axes of our Milky Way galaxy 
model. The Galactic plane corresponds to the XY projec- 
tion. Distances are in kpc. 



flipped with respect to the disc plane. Two of the orbits 
for the A/'-body satellites are shown in Fig. |2| All stars 
of a given dwarf galaxy are assumed to be of the same 
age. Masses are drawn from a mass-function and My 
and { V — I) are derived fr om a low-metallicity isochrone 
rfrom lGirardi et al.Ll200Ql) of the appropriate age. 



4. MODELLING THE GAIA SURVEY 



The 6 phase space coordinates (r, v) for each star in the 
Galaxy and satellite model are referred to the Solar po- 
sition and velocity and converted to the 5 astrometric 
parameters; position (^,6), parallax tu, proper motions 
M£* = M£ cos h and /i^, and the radial velocity Vrad- The 
conversion is done using the standard prescriptions for 
transforming Cartesian position and velocity coordinates 
into astrometric parameters and radial velocities, such as 
described in in Volume 1, Section 1.5.6, of the Hippar- 
cos Catalogue (ESA, 1997). The astrometric errors are 
added as a function of G (Gaia broad band magnitude) 
and {V -I): 
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105z+ 1.3^2 + 6 X 10 
x[0.96 + 0.04(V-/)] 

where z = iQ^-'^iG-i^) 
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G= V + 0.51 - 0.50 X ^06 + 
0.065 x{V-I- 
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The parallax errors for the simulated Gaia data of the 
Milky Way model are shown in Fig. |3l The mean po- 
sition and proper motion errors are ctq = 0.87cr^ and 
c^/i = 0.75cr^, respectively, and the variation of the as- 
trometric errors with ecliptic latitude (3 is included. The 
radial velocity errors for OB A- type stars are 0.25, 4, 10, 
50 km s~^ at V = 10, 15, 16, 17, respectively. For 
FGKM-type stars the errors are 0.1, 1, 2, 6 km s~^ at 
V = 10, 16, 17, 18. For stars fainter than V = 17 (OBA) 
or V = 18 (FGKM) radial velocities are not available. 
These error prescriptions are according to the Gaia Con- 
cept and Technology Study Report lESA>,2000t) . 
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Figure 3. Parallax errors vs. V -magnitude in the simu- 
lated Gaia catalogue for 1 million stars from the model 
Milky Way. The grey scale indicates the number of stars. 
The spread in the values at each V is caused by the colour 
dependence of the errors (see Eq.^). 



4.1. Combining the Simulated Milky Way and Satel- 
lite data 

Having made a considerable effort to realistically sim- 
ulate the number of Galaxy stars that is expected to be 
seen by Gaia, we want to ensure that the satellite simula- 
tions are properly added to the Galaxy data. This means 
that the number of satellite particles in our simulated Gaia 
catalogue should be a realistic fraction of the number of 
Galactic particles. Getting this right is not trivial and we 
explain our solution to this problem here. 

Given a certain distribution of stars along the orbit of 
a particular dwarf galaxy (corresponding to an A^-body 
'snap-shot'), the number of stars from this satellite that 
will end up in the Gaia catalogue depends on three fac- 
tors: 1. The overall number stars (i.e., luminous parti- 
cles) in the dwarf galaxy. This number is determined by 
its overall luminosity and the stellar mass function. 2. 
The Gaia survey limit leads to an upper limit My^max on 
the absolute magnitude of visible satellite stars. 3. The 
variation of My^max along the satellite orbit, caused by a 
variation in distance from the Sun. 

At each distance Si where an N-bo&y particle is located 
the fraction of visible satellite stars fi can be calculated: 



fi{Si 



m(Mv, max) 



(3) 



where Mv,maA^i) Mim - Slogs^ + 5, while ^(m) 
is the mass function and m^p its upper limit. The 
overall fraction / of visible A/'-body stars is / = 
J2i fii^i) /^sim- As Fig.|4|illustrates the overall fraction 
of visible stars for a fully populated mass function can be 
very small depending on the distribution of satellite stars 
along the orbit. In addition we are faced with the prob- 
lem that real dwarf galaxies contain up to 10^ stars but 



Figure 4. fi vs Si for satellite model 1. The shaded his- 
togram (right vertical scale) shows the distribution ofsi. 
The solid line (left vertical scale) shows fi{si). Features 
in this line reflect features in the luminosity function. The 
dashed line shows the fi{si) when assuming all N-body 
stars to be brighter than the faintest satellite star that en- 
ters the Gaia survey. We calculated fi using a 10 Gyr 
isochrone with Z = O.OOA and (_{m) ex m~^-^. The over- 
all fraction of visible N-body particles is only 2.5 x 10~^ 
(I) and 10~'^ for the two cases considered. 



our A'-body models contain 1 million particles only. Of 
these 'expensive' A/'-body particles we want to waste as 
little as possible. An obvious step is to assume that all 
N-body particles represent stars that are brighter than the 
faintest star that can enter the Gaia survey given the dis- 
tance distribution of the A^-body particles. This will raise 
the overall visible fraction A/'-body particles as illustrated 
by the dashed line in Fig.|4| 

This hints at the following solution: assume that all A^- 
body particles represent a bright tracer population, such 
as AGB stars. That is. My < My^tracer for all par- 
ticles. This will further raise the overall visible frac- 
tion /. The details of thi s procedu re are discussed in 
Brown. Velazquez. & Aguitol (l2QQ4l) where it is shown 
that depending on the total luminosity (mass) of the sim- 
ulated satellite one can then easily retain the majority 
of the A/'-body particles while still preserving the varia- 
tion of the visible fraction of tracer stars along the debris 
stream. 



TRACING DEBRIS STREAMS IN THE GAIA 
CATALOGUE 



Various methods have been proposed and used to recover 
satellite remnants from surveys of Galactic phase space. 
Figure|5lshows an example of energy vs angular momen- 
tum (E-Lz) diagrams. For this figure we combined 18 
million stars from the Monte Carlo model of the Galaxy 
with the satellite models 1 (at 10 Gyr) and 4 (at 5 Gyr). 
The value of My^tracer is 0.4 and 1.0 respectively. Keep 
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Figure 5. Energy vs angular momentum diagrams. The top row shows the true (simulated) E-L^ distribution and the 
bottom row the energy and angular momentum as they would would appear when derived from the Gaia catalogue. Erom 
left to right the panels show the data for the Monte Carlo model of the Milky, for the two satellites (models 1 and 4 at 10 
and 5 Gyr, respectively), and the combined data. Only stars with positive parallaxes and available radial velocities are 
shown. Note how the compact E-L^ distribution of the debris streams turns into a more complex structure in the lower 
diagrams due to the way the astrometric and radial velocity errors propagate. The parallax errors (see Eig.^ are by far 
the dominant contribution to these structures. 



in mind that the figures under-represent the real contrast 
between Galaxy and satellite. 



6. FUTURE WORK 



We intend to use this work to study in detail the retrieval 
of debris streams from the Gaia catalogue by making use 
of the accurate phase space data that Gaia will provide 
throughout the Galaxy. However, as can be appreciated 
from the E-Lz diagrams shown in Fig. 13 it will not be 
possible to trace complete debris streams based on phase 
space data alone. Fortunately, using the photometric in- 
formation from Gaia will provide astrophysical parame- 
ters for the different stellar populations. This will enable 
a drastic narrowing down of the parts of the Gaia cat- 
alogue that need to be searched. Investigating this will 
require improvements to the simulations. 

The simulation of actual Gaia photometry from the broad 
and medium band photometers is possible using the tools 
developed at the University of Barcelona (see the contri- 
butions by Jordi et al. and Carrasco et al. in this volume). 
An extinction model can be included in an approximate 
way using currently available 3D extinction models such 
as the one presented by Drimmel et al. (this volume). Fi- 



nally, if in reality the Galactic halo has not had time to 
relax and wipe out the clumpy remnants of past merger 
events, it will be necessary to replace the smooth halo 
component in our Milky Way model with a clumpy one. 
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